Counting statistics of photons produced by electronic shot noise 
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A theory is presented for the photodetection statistics of radiation produced by current fluctu- 
ations in a phase-coherent conductor. Deviations are found from the Poisson statistics that would 
result from a classical current. For detection in a narrow frequency interval 6u), the photocount 
distribution has the negative-binomial form of black-body radiation if eSui is less than the mean 
current I in the conductor. When electronic localization sets in, I drops below eSuJ and a different 
type of super-Poissonian photon statistics results. 
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Some interesting recent developments in mesoscopic 
physics have arisen from the interplay with quantum op- 
tics Q. To mention two examples, the Hanbury-Brown- 
Twiss effect for photons is inspiring the search for its elec- 
tronic counterpart B , while single-electron tunneling has 
been used to create a single-photon turnstile device ]|). 
An appealing subject of research, in line with these devel- 
opments, is the study of a mesoscopic conductor through 
the quantum optical properties of the radiation produced 
by the current fluctuations. It is a textbook result [Q, 
due to Glauber || , that a classical current produces pho- 
tons with Poisson statistics. What is the photon statis- 
tics for a fully phase-coherent conductor? That is the 
fundamental question addressed in this paper. 

It is a timely question in view of a recent proposal 
by Aguado and Kouwenhoven Q to use photon-assisted 
tunneling in a device containing two quantum dots in 
series as a detector for the microwave radiation emitted 
by a nearby mesoscopic conductor (e.g. a quantum point 
contact) . One such device by itself can only give informa- 
tion on the mean rate of photon production, calculated 
in Refs. HQ], but a pair of devices could measure the 
time-dependent correlations and hence could detect de- 
viations from Poisson statistics due to photon bunching 

We will calculate these fluctuations for an idealized 
model of a photodetector, the same model that leads to 
the Glauber formula of photodetection theory || . In this 
formula the photocount distribution is expressed as an 
expectation value of normally ordered photon creation 
and annihilation operators. (Normal ordering means that 
all creation operators are brought to the left of the annihi- 
lation operators.) We will see that the ordering inherited 
by the electron current operators involves not only a nor- 
mal ordering, but in addition an ordering of the incoming 
current with respect to the outgoing current. 

We present a general formula for the variance of the 
photocount in terms of the transmission and reflection 
matrices of the conductor. A particularly simple result 
is obtained in the limit that the frequency interval 8uo 
of the detected radiation is small compared to the mean 
(particle) current I/e through the conductor: The photo- 
count distribution P(n) for a long counting time t is then 



proportional to the negative-binomial coefficient ( n+ ^ _ ) 
(with v — tScj/2tt ^> I) pcj . This is the photonic coun- 
terpart of the (positive) binomial counting distribution 
for electrical charge Jll|-Q. In the localized regime the 
condition Slu I/e breaks down and a different non- 
Poissonian distribution results. 

Starting point of our analysis is an expression for the 
photocount distribution as a time-ordered expectation 
value of the electric field operator M , 



P(n) = (T ± -W n e- w ), 



(1) 
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Here P(n) is the probability to detect n photons in a time 
interval r and E{t) is the operator of the detected mode 
of the electric field in the Heisenberg picture. (We as- 
sume for simplicity that a single mode is detected.) The 
detector has sensitivity a(w) at frequency to > 0. The 
symbol T± indicates the Keldysh time-ordering of the 
time-dependent operators: times i_ to the left of times 
t+, earlier i_ to the left of later t_, earlier t + to the right 
of later t + . 

The Glauber formula mentioned in the introduction is 
obtained from Eq. (|l|) by substituting the free-field ex- 
pression 

poo 

E llcc (t)(x du (a) (u)e luJt + a{u)e- lu)t ) (2) 
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for E(t) and by making the rotating-wave approxima- 
tion [neglecting e^" +a) )*, retaining e^ w_ " )*]. The time 
ordering of the electric field then becomes normal order- 
ing of the photon operators ad , a. Our goal is instead to 
go from Eq. (Q) to an expression in terms of the electron 
operators c\c that constitute the current operator I(t). 

The electron and photon degrees of freedom are cou- 
pled in the Hamiltonian via a term — J dr j(r, t) ■ A(r, t), 
where j is the electron current density operator and A 
the electromagnetic vector potential. This bilinear cou- 
pling leads to a linear integral relation between E and 
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E(t) = E ilee (t) + / dt'g(t - t')I{t'). 



(3) 



The propagator g{t) vanishes for t < because of causal- 
ity. (We are neglecting retardation of the electromag- 
netic radiation.) We assume that electrons and photons 
are uncoupled for t — * — oo, the photons starting out in 
the vacuum state. Substitution of Eq. (||) into Eq. (|l|) 
leads to a correlator involving the non-commuting oper- 
ators Ef TCC and I. Flcischhauer has shown jl5| that the 
vacuum term E{ ree may be removed from the correlator 
if the Keldysh time-ordering of E is carried over to I: 



W = 



duja{uj) // dfdf'e^"-^ 



dt-dtA 



xg(t' -t_)g(t" -t + )I{t_)I{t + ). 



(4) 



To leading order in g we may neglect the coupling to 
the photons in the time dependence of I(t). For free 
electrons the time dependence is given by 16 
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dede' e i(£ - £ ' V( £ )M(£, e')c(e'), (5) 



M(e,e') = SHe)DS{e') - D, 



D = 




1 



s = 



r' t' 



(6) 
(7) 



using units such that Ti = 1 and e = 1. We have intro- 
duced the scattering matrix S(e) (with N x N reflection 
and transmission submatrices r,r',t,t') of the N propa- 
gating modes at energy e (relative to the Fermi energy at 
e = 0). The scattering geometry is illustrated in Fig. ||. 
The detection matrix D selects the current in one of the 
two leads, arbitrarily chosen to be the right lead in Eq. 
(0) . The total current I = 7 out — Ii n is then the difference 
of the current / ut coming from the left and the current 
/in coming from the right. These two currents Tout and 
Ii n are defined as in Eq. (|^) , with the matrix M replaced 
by S^DS and D, respectively. 




FIG. 1. Illustration of the scattering geometry studied in 
the text. An electrical current flowing through a constric- 
tion emits microwave radiation that is absorbed by a nearby 
detector. 



The separation of / into outgoing and incoming current 
operators is convenient because they have simple commu- 
tation relations: (i) I ou t{t) commutes with I ut(t'); (ii) 
ii n (t) commutes with L m (t'); (hi) I ou t{t) commutes with 
iin(i') if t < t'. It follows that Keldysh time-ordering of 
the current operators is the same as an ordering whereby 
the operators /i n (i_) are moved to the left and I ln (t+) 
to the right of all other operators — irrespective of the 
values of the time arguments. 

Now that we have liberated ourselves from the time 
ordering we are free to take Fourier transforms, 

1 r°° 

P(n) = (0-W n e- w ), W= duja(uj)U\uj)U(cu), 



U(u) 



00 du' 
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K(lu - uo')g{J)[I ont (Lo') - Z ln (u/)]. (8) 



The Fourier transforms of g,/ utj^in & re defined as 
f{u>) = dt e lujt f (t) , and we have abbreviated 

K{lo) = JJ" dt' e luJt . The symbol O indicates ordering 
of the current operators from left to right in the order 
J L 4ut > 4ut. An- According to Eq. (|), the current oper- 
ators are related to the electron creation and annihilation 
operators by 



7 ou tM = j dec\e)S\e)DS(e + uj)c(e + u>), (9a) 
I in (u,) = j deJ{e)Dc{£ + w), (9b) 



Eqs. (|^) and (|j) form the required relation between the 
photocount distribution and the electron creation and 
annihilation operators. 

The mean photocount n has been studied before, by 
Lesovik and Loosen [fjj and by Aguado and Kouwenhoven 
To make contact with that work we take the exper- 
imentally relevant limit of a long detection time r. We 
may then discretize the frequencies as oj n — n x 2tt/t. 
In this discretization the kernel K becomes a Kronecker 
delta, K{u n - Lo m ) = r5 nm , hence U(uj n ) = g(uj n )I(uo n ). 
The factorial moments 



{n p ) t = n(n - l)(n - 2) • • • (n - p + 1) = {OW p ) (10) 

of the distribution (||) in the long-time limit take the form 
lE7| 



K) f = (0[/do;7M/ t MlMF), 



(11) 



with ~f(uj) = a(ui)\g(uj)\ 2 . For the first moment the or- 
dering operator O can be omitted and we find 



oo 

// = / dw7(w)(/ t (w)/(u))), 
'o 
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(Jt («;)/(«;)> =t dte tut (I(0)I(t)), 



(12) 
(13) 



in agreement with Refs. 
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For the double-quantum-dot photodetector of Aguado 
and Kouwenhoven Q the response function j(oj) is 
sharply peaked at the frequency 17 of the inelastic tunnel- 
ing transition, with a width 5u> <C f2. Its integrated mag- 
nitude j(n)8uj ~ (ZT /Q) 2 depends on the impedance Z 
(in units of h/e 2 ) of the inductive coupling and the ratio 
of f2 and the tunnel rate T between the quantum dots. 
(Typically, Z < 1 and V < Q, so that jSoj < 1.) In 
that device the inelastic transition can be either upwards 
or downwards in energy, corresponding, respectively, to 
the absorption or emission of a photon. Here we consider 
only the case of photodetection by absorption, which is 
the relevant case for the study of shot noise in the con- 
ductor @. 

We now go beyond the first moment to study the en- 
tire photocount distribution. We note that P(n) in Eq. 
(||) would be simply a Poisson distribution, 
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Poisson (^) — 71 e , 

n! 



if the current / would be a classical quantity instead of 
a quantum mechanical operator. This is in accordance 
with Glauber's finding Q that the radiation produced by 
a classical current is in a coherent state (since a coherent 
state has Poisson statistics) . 

To find the deviations from the Poisson distribution 
due to quantum statistics, let us consider the case of a 
conductor connecting two electron reservoirs in thermal 
equilibrium at temperature T. The system is brought 
out of equilibrium by application of a voltage difference 
V between the left and right reservoirs. Expectation val- 
ues are given by the Fermi function /(e) = (e £ / T + 



(4(e)c j (e'))=6 ij 5(e-e')f(e-fa), 



(15) 



with higher order expectation values obtained by pair- 
wise averaging. The potential fa equals V for the left 
reservoir (mode indices i — 1,2,..., N) and for the 
right reservoir (i = N + 1, N + 2, . . . , 2N). 

For simplicity we restrict ourselves to zero tempera- 
ture, when f{e — fa) becomes the step function 9(fa — e) 
(equal to 1 for e < fa and for e > fa). The mean and 
variance of the photocount arc then given by 



=- / dujV, 

ZTV 



(16) 
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Varn = nH / dujSf 2 H // duduj' 7(^)7(0/) 
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(17) 
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Jo 

JV = -f(iv) / deTiT e p £ - LU , 

■J UJ 



(18) 
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x9(e-Lu-uj'), (19) 

A 2 = TePe-uTePe-wOie ~ Cj)0( £ - W 7 ), (20) 

A 3 = p e T e+u p e T s+u ,6(V - e - u)6(V -e- uj'). (21) 



We have abbreviated r £ = t(e)t'(e) — 1 — p e . 

The formula (16h for the mean photocount is known 
the result ( |17| ) for the variance is new. The first 
term n on the right-hand-side corresponds to Poisson 
statistics. The other terms describe the excess noise, 
consisting of one term containing a single integral over 
frequency and three more terms containing double fre- 
quency integrals. For narrow-band detection the single 
frequency integral dominates. More precisely, if 7(0;) is 
non-zero in a narrow frequency range 8u> <C V, then 



Varn = (t5uj /2n)M(l + N). 



(22) 



The correction terms from the double frequency integrals 
are smaller by a factor 8uj / 1, with / oc V Tr W the mean 
electrical current flowing between the reservoirs. 

In this regime of narrow-band detection one can also 
calculate easily the higher order moments of the photo- 



(14) count. The factorial cumulants ((n p ))f are given by 



(( n P)) { = {T5u}/2ir)(p-l)\Ar p . 



(23) 



The probability distribution P(n) can be reconstructed 
from the factorial cumulants via the generating function 
F{0 = E^Li(£ p /p!)((™ p ))f, by means 01 the formula 



1 r1 n 

P(n) = -, lim ^-e F ^. 



(24) 



The probability distribution corresponding to Eq. (|23|) is 



P(n) 



n + v - 1\ 

n ) (l+Af) n +»' 



(25) 



which is the negative-binomial distribution with v = 
tSuj/2-k degrees of freedom. (For non-integer v the bi- 
nomial coefficient should be interpreted as a ratio of 
Gamma functions.) It approaches the Poisson distribu- 
tion ( fl4"| ) in the limit M — > 0, v — > 00, at fixed n = vN . 
The negative-binomial distribution is known in quantum 
optics as the distribution of black-body radiation jHj . The 
role of M is then played by the Bose function [e w ' — 1) _1 
at the temperature of the black body. In both contexts M 
is a small parameter, and hence the corrections to Poisson 
statistics are small: for a black body, Af is small because 
uj ^> T, while for the electrical conductor J\f ~ 7J <C 1. 

We have seen that the negative-binomial distribution 
results if contributions of order Soj/I can be neglected. 
If the electrical conductor is metallic, it is sufficient that 
5uj V, since the conductance (in units of e 2 /h) is 
greater than 1 in a metal. In the localized regime, on the 
contrary, the conductance becomes exponentially small 
and terms of order Suj/I start playing a role — even if 
Soj <C V. To illustrate this, let us assume that the trans- 
mission through the conductor is so small that only terms 
linear in t £ need to be retained. This leaves only the term 
At in Eq. (Il7|), so that 
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Varn = n+ — (jSuj) 2 6(V - 2w) / cfeTrr E . (26) 

More generally, the factorial cumulants are given by 

((n p h = ^{i5loY9{V - pu) [ de Trr £ . (27) 

The full distribution P(n) can be reconstructed by means 
of the inversion formula (p^), but does not have a sim- 
ple closed-form expression. We note that the deviations 
from Poisson statistics are again small because jStu <C 1 . 

Much larger deviations can be obtained if coherent ra- 
diation from a reference source at frequency Q is superim- 
posed prior to detection. Such homodyne detection not 
only amplifies the deviations from Poisson statistics, it 
also provides a way to measure the counting distribution 
of electrical charge. To see this, we note that homodyn- 
ing amounts to the replacement of the current operator 
I(oj) by J (a;) + IqS(uj — ft), where Iq is some known clas- 
sical current (a c-number, not an operator). For Iq 3> I 
we find from Eq. ( |Tl| ) for the factorial cumulants of the 
photocount distribution the expression 

K)) f = Mr 7 / 2 /2^) + (27/0)^)), (28) 
where ((Q p )) is defined through the generating function 

Comparison with Refs. pT| , p^ | shows that {(Q p )) has the 
interpretation of the cumulant of the charge Q transmit- 
ted through the conductor (in units of e). Levitov and 
Lesovik Q| proposed to measure the charge counting 
distribution from the precession of a spin i coupled to 
the current. The photodetection scheme proposed here 
provides an alternative, and possibly more practical, way 
to count the charge without breaking the circuit. 

In summary, we have presented a solution to the classic 
problem of the statistics of radiation produced by a fluc- 
tuating current. We go beyond the textbook result by 
considering a fully phase-coherent conductor, and find 
small deviations from the Poisson statistics associated 
with a classical current source. The deviations might 
be measured using an array of double-quantum-dot pho- 
todetectors Q. The deviations can be amplified by ho- 
modyning, in which case they are directly related to the 
statistics of the electrical charge transmitted through the 
conductor JTTJ . We have given specific results for a con- 
ductor between normal reservoirs in thermal equilibrium, 
but our general formulas can be applied to more special 
current sources as well. The applications to entangled 
Jl9| or superconducting electrons seem particu- 

larly interesting. 
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